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Abstract 

Molecular nanomagnets, besides promising to open new frontiers in technology, have attracted 
huge interest in the scientific community because they can exhibit the phenomenon known as 
quantum tunnelling of the magnetization, i.e. coherent fluctuations of the direction of the total 
spin vector. In this paper we study a different quantum phenomenon involving fluctuations of the 
magnitude of the total spin vector. These fluctuations are related to the mixing between states 
with different spin quantum number, and imply new macroscopic effects, which we theoretically 
investigated in the Mn-[3x3] grid. 
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Molecular nanomagnets (MNMs)^^ are clusters containing a finite number of transition- 
metal ions whose magnetic moments (spins) are so strongly coupled that at low tempera- 
ture each molecule behaves like a single-domain particle with fixed total spin. Being at the 
crossover between classical and quantum regimes, MNMs exhibit at the same time classical 
properties of macroscopic magnets such as magnetization hysteresis, and quantum phenom- 
ena like tunneling of the direction of the total spin through energy barrier a 4 ! 5 ' 6 . MNM 
systems are interesting also for potential technological applications, as envisaged for the im- 
plementation of quantum computing algorithms 7 , or for dense and highly efficient memory 
devices^. Here we study a new macroscopic manifestation of a quantum phenomenon involv- 
ing fluctuations not only of the direction, but also of the magnitude of the total spin, and we 
show that it is realized in a Mn-[3x3] grid^. Recognizing the effects of these fluctuations is 
essential to achieve a satisfactory understanding of the role played by quantum mechanics 
in complex mesoscopic magnetic systems. Quantum magnetic phenomena were identified 
even in molecules of great biological interest such as ferriti n 10 ! 11 . 

The advantage of studying quantum phenomena in MNMs is that the vanishingly small 
interaction between different molecules allows single-molecule phenomena to be observed at 
a macroscopic scale, because the crystal behaves like a collection of independent objects^ 1 ^, 
each described in general by spin Hamiltonians of the form^ 



where Sj are spin operators of the i th ion in the molecule. The first term is the isotropic 
Heisenberg exchange interaction, the second and third terms describe the local crystal- 
field and the anisotropic intra-cluster spin-spin interactions. The last term is the Zeeman 
coupling with an external field B in which isotropic ^-factors are assumed. While the 
Heisenberg term is rotationally invariant and therefore conserves the length |S| of the total 
spin S = Y2i s i> the anisotropic terms break rotational invariance and do not conserve this 
observable. Nevertheless, since the Heisenberg contribution is usually largely dominant, |S| 
is nearly conserved, and the energy spectrum of H consists of a series of level multiplets with 
an almost definite value of | S | (expressed in terms of the quantum number S as a/ S(S + 1)) . 
Thus, quantum fluctuations of |S|, which are associated with mixing of states with different 
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value of the quantum number S (" S- mixing"— ) , either are zero or are expected to produce 
negligible effects on the macroscopic behavior, and are therefore neglected in virtually all 
studies. A major theoretical goal would be to identify a clear macroscopic signature of such 
fluctuations, and a model system displaying such effect. 

In this paper we show that favorable conditions are met in the Mn(II)-[3x3] grid-like 
cluster [Mn 9 (2POAP-2H) 6 ](C10 4 ) 6 -3.57MeCN-H 2 (hereafter Mn-[3x3], see Fig. 1), a sys- 
tem that has been recently characterized by magnetization and torque measurements*^. 
Our theoretical calculations suggest grid-shaped molecules as good candidates to study fluc- 
tuations of |S|, since the lowest level is expected to display significant S-mixing with the 
first excited multiplet. In fact, unlike in ideal ring-shaped molecules^, in the [3x3] grid 
the two lowest manifolds of the exchange part of the spin Hamiltonian belong to the same 
irreducible representation of the molecular symmetry group. Therefore, the application of a 
suitably oriented magnetic field induces a series of anticrossings (ACs) between the ground 
state and levels originating from higher excited manifolds (see red lines in Fig. 2). At the 
AC fields, ^-mixing in the ground state is maximum, and quantum fluctuations of |S| are 
greatly enhanced. Torque^ and neutron^ 7 - experiments on Mn-[3x3] show that the zero-field 
gap between the two lowest S-multiplets is small enough for the ACs to occur at fields within 
experimental reach. 

Mn-[3x3] crystallizes in the space group C2/C, and the cation [Mn 9 (2POAP-2H) 6 ] 6+ ex- 
hibits a slightly distorted S4 molecular symmetry with the C2 axis perpendicular to the plane 
of the cluster^. The average distance between the Mn(II) ions is 3.93 A, the smallest distance 
between clusters is larger than 8 A. For a cluster composed of nine interacting Mn(II) spins 
with Sj = 5/2 the dimension of the Hilbert space is 10077696. The difficulties related with 
this huge dimension have been overcome by exploiting both the irreducible tensor operator 
technique and the Lanczos algorithm for the exact diagonalisation. The two-step procedure 
already developed 1 -^ has allowed the inclusion of S'-mixing effects in the calculation. The 
exchange integrals and the single-site and spin-spin anisotropy tensors have been determined 
by neutron spectroscopy 17 . All exchange integrals between nearest-neighbors are found to 
be nearly equal to 0.47 meV, apart from J 18 , J 78 , J 34 and J 45 (see Fig. 1), whose values are 
0.33 meV. Next-nearest-neighbor exchange interactions can be neglected. Concerning the 
local crystal-field, the second term in Eq. 1 can be approximately rewritten as 
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Si ■ Di • * = D [4 " + !)]> (2) 

i i 

with D = —6.1 /xeV. The intracluster dipole-dipole interaction Djj has been evaluated 
within the point-dipole approximation^. At last, gi = 2 is assumed, as appropriate for 
Mn(II) ions. 

With these experimentally determined parameters, the energy spectrum for B = con- 
sists of many level multiplets with an almost definite value of S. The ground multiplet has 
S = 5/2, and the four lowest-lying excited multiplets have (in order of increasing energy) 
S = 7/2, S = 3/2, S = 3/2, S = 9/2. These multiplets are separated by the isotropic 
exchange and split by the anisotropic interactions. The lowest levels are shown in Fig. 2. 
The application of a magnetic field B in a direction outside the grid plane and different from 
that of the axis produces several ACs involving levels belonging to different multiplets 
(see Fig. 2). As the AC fields B c are approached the multiplet mixing is enhanced. For 
B = B c , the spins in each cluster oscillate coherently between states with different values of 
S, which therefore is no longer a good quantum number. This can be inferred for example 
for the ground state ACs (indicated by the arrows in Fig. 2) from the inset in Fig. 2, which 
shows the field-dependence at T = of the quantity S e ff, defined through the relation 
S e ff(S e ff + 1) = (S 2 ). When B is close to B c the value of S e ff is intermediate between 
half integers, e.g. 5/2 and 7/2 at the first level AC, thus confirming that the ground-state 
wavefunction is a superposition of different total-spin states. 

This opens a scenario in which not only the direction of the total spin fluctuates in time, 
as in the case of the well-studied quantum tunnelling of the magnetization^^, but even 
its length fluctuates (quantum dynamics of the total spin, or briefly QDTS). Moreover, by 
properly tuning the direction of the applied magnetic field, the AC splitting can be made 
as large as several Kelvins, thus overcoming the problem of decoherence due to hyperfine 
fields and to cluster-cluster interactions. Most importantly, in the Mn-[3x3] grid several of 
these ACs involve the two lowest energy levels, thus opening the possibility to detect the 
QDTS by means of macroscopic low-temperature magnetic bulk techniques. In particular, 
our calculations show that at low temperature, in correspondence to each AC between the 
two lowest lying states, a sharp peak should be detected in the torque signal (as function 
of the field intensity B) when the direction of the field 9 is close to the grid plane (8 = 0°) 
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or almost perpendicular to it {6 = 90°) (see Figs. 3 and 4). The precipitous drop of the 
torque signal as a function of 9 near 8 = 0° and 9 = 90° is due to a change of sign imposed 
by symmetry (as M x B = when B points along a principal direction of the susceptibility 
tensor of the system). 

By comparing Figs. 3 and 4 with the corresponding spectrum reported in Fig. 2, it 
is evident that a torque peak appears in correspondence to each AC involving the ground 
state. While the first peak at low field corresponds to an intra-multiplet AC, the peaks at 
about 6.6, 8.9 and 11.2 Tesla represent a direct consequence of the coherent superposition 
of different total spin quantum states. The link between the peaks in the torque signal 
and the level ACs is demonstrated in Fig. 4, where the torque curves calculated with the 
mixing artificially forced to zero are shown by colored dashed lines. When S-mixing is 
neglected, every inter-multiplet AC becomes a crossing, the fluctuations of the total spin 
vector are suppressed and the peaks in the torque curve disappear. Thus, without S'-mixing 
the torque curve should exhibit the step-like field dependence usually observed in MNMs^. 
The comparison between solid and dashed curves in Fig. 4 shows clearly that the predicted 
peaks in the torque would be a direct macroscopic manifestation of the QDTS at the level 
ACs. Also, these spin fluctuations and the associated peaks are of purely quantum origin, 
as these are absent in the classical version of the Hamiltonian Eq. (1) (see the gray dashed 
curve in Fig. 4). 

The fact that these torque peaks directly reflect the large enhancement of quantum fluc- 
tuations of | S | at the ACs can be understood by a simple physical picture. In the following 
we define the z' axis parallel to B. The torque is proportional to the magnetic response 
perpendicular to the field direction, i.e. to (S x >). Fluctuations of |S| are accompanied by 
fluctuations of the magnetization (S z >) 2 ®, and these latter are deeply connected with (S x '). 
Accordingly, S x > tracks the increase and decrease of these fluctuations while sweeping over 
the AC, leading to a peak in the torque. Indeed, near the ACs, the ground-state wavefunc- 
tion can be written to a good degree of approximation as 

\G) = -±=(a{B)\r 1 ,M) + b{B)\r 2 ,M+l}), (3) 

where a(B) 2 + b(B) 2 = 2 by normalization. For B < B c , b(B) ~ 0, a(B) ~ y/2, and 
the ground state reduces approximately to an eigenstate of S Z ', |Fi,M), where M is the 
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corresponding eigenvalue and Ti represents the set of additional labels necessary to identify 
the state. For B ^> B c , a(B) ~ 0, b(B) ~ ^/2, and the ground state reduces approximately 

to |r 2 ,M + i). 

At T = the torque is given by 

t oc B(S X >) = Ba(B)b(B)(T lj M\S x ,\T 2 ,M +1). (4) 
Since (Ti, M\S X >\T 2 , M + 1) is almost independent on B, the field dependence of r comes 



entirely from the factor Ba(B)b(B) = Bb(B)^/2 - b(B) 2 . 

Quantum fluctuations of the magnetization are usually characterized by the quantity 

(AS Z >) = {S z ,} - {S z >) = — — . (5) 

Therefore, 

r oc 2BAS Z /. (6) 

In case of anticrossing, AS z r is maximum at B c , where b(B c ) = a(B c ) = 1, and the torque 
peaks. In case of crossing AS Z > is always zero because either b{B) = or a(B) = 0, and the 
torque does not peak. As a further confirmation of this picture we show in Fig. 5 numerical 
calculations of the fluctuations of S z > ("parallel") and of S x > (" perpendicular"), with and 
without S — mixing. 

The torque is the macroscopic magnetic bulk technique displaying the clearest signature 
of S'-mixing. When the temperature T increases these effects wash out, and indeed published 
torque measurements on Mn-[3x3} 9 do not display peaks because data were collected at too 
high a temperature. 

We have also calculated the low-T heat capacity as a function of the applied field (see 
Fig. 6). This quantity provides a way to assess the value of the AC g function of 

the field direction. For 6^0 anticrossings open up and the heat capacity does not vanish 
at B c . 

In this paper we have identified a new clear-cut macroscopic manifestation of a quantum 
phenomenon characterizing the microscopic dynamics of magnetic clusters. We have shown 
that torque measurements provide direct evidence of quantum fluctuations of the total spin 
length, and we have studied a model system displaying such effect. The comprehensive 
understanding of the spin dynamics in mesoscopic systems like molecular nanomagnets is 
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essential to understand the crossover between quantum and classical mechanics and to make 
rational design of these compounds possible, especially in view of the promising applications. 
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Fig. 1: A schematic representation of the molecular structure of the Mn-[3x3] grid 
showing the network of manganese (blue) and oxygen (green) bonds. Other atoms are 
omitted for clarity. The red arrows indicate the spin directions in the classical ground state. 



Fig. 2: Calculated energy levels of the Hamiltonian Eq.(l) with the parameters given in 
the text, and with the direction of the applied field forming an angle 9 = 2.8° with the grid 
plane. Energies are plotted as functions of the applied field intensity B. The ground state 
energy is set equal to zero. Arrows indicate the anticrossings produced by S'-mixing between 
the two lowest levels (highlighted in red). The inset shows S e ff as functions of B for the 
same angle at T = 0, where S e ff is defined through the relation S e ff(S e ff + 1) = (S 2 ). 



Fig. 3: Two-dimensional plot of the calculated torque as a function of the applied field 
direction 9 and intensity B at T = 0.05 K with parameters of Hamiltonian Eq.(l) as given 
in the text. 



Fig. 4: Plots of the calculated torque at T = 0.4 K as a function of the applied 
magnetic field B for two selected directions 9 lying close to the grid plane. Colored dashed 
lines represent the calculated torque at the same angles with S'-mixing eliminated by hand, 
i.e. with quantum fluctuations of |S| removed. The dashed gray line represents the torque 
calculated using the classical version of the Hamiltonian Eq. (1) at T = for 9 = 2.8°. The 
classical value of the torque has been rescaled (by about 3) to fit the figure. The inset shows 
the calculated torque at T = 0.4 K as a function of the applied field intensity B for several 
directions 9. 



Fig. 5: Numerical calculations of the T = fluctuations of S z > ("parallel") and of S x > ( 
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perpendicular"), as a function of the applied magnetic field (forming an angle 9 = 2.8° with 
the grid plane), with and without S — mixing. Parameters are given in the text. 



Fig. 6: Numerical calculations of the heat capacity as a function of the applied magnetic 
field for several directions 9, at T = 0.4 K. 
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